Adsorption in non interconnected pores open at one or at both ends: 
a reconsideration of the origin of the hysteresis phenomenon. 
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We report on an experimental study of adsorption isotherm of nitrogen onto porous silicon with 
non interconnected pores open at one or at both ends in order to check for the first time the old 
(1938) but always current idea based on Cohan's description which suggests that the adsorption of 
gas should occur reversibly in the first case and irreversibly in the second one. Hysteresis loops, the 
shape of which is usually associated to interconnections in porous media, are observed whether the 
pores are open at one or at both ends in contradiction with Cohan's model. 



PACS numbers: 68.43.-h, 61.43. Gt, 64.70.Fx, 05.70.Np 



Adsorption isotherms in mesoporous materials exhibit 
two main typical features: (i) a sharp increase in the 
amount of gas adsorbed at pressure below the saturated 
vapor pressure of the bulk gas which is attributed to cap- 
illary condensation in pores of confining sizes; (ii) the 
desorption process fails to retrace the path of the ad- 
sorption one, so that a hysteresis loop appears revealing 
the irreversibility of the phenomenon. 

The first explanation of the capillary condensation was 
proposed by Zsigmondy JlJ who, for this purpose, refered 
to the macroscopic Kelvin equation. In order to explain 
the hysteresis phenomenon, Cohan ||] went further and 
proposed different scenari for filling and emptying of a 
pore which he also described with the macroscopic con- 
cept of meniscus. In a cylindrical pore open at both 
ends, the shape of gas/adsorbate interface is different 
during the adsorption and desorption processes leading 
to a hysteresis phenomenon. On the contrary, in a cylin- 
drical pore closed at one end, the adsorption is reversible 
since the same meniscus is present during the filling and 
the emptying of the pore. 

The irreversibility of the adsorption process in a cylin- 
drical pore open at both ends, expected by Cohan on 
the basis of thermodynamic considerations, has been 
predicted by new numerical approaches. Thus, Grand 
Canonical Monte-Carlo (GCMC) simulations analyt- 
ical approaches Q] together with Density Functional The- 
ory (DFT) H have also proposed that the hysteresis phe- 
nomenon is an intrinsic property of the phase transition 
in a single idealized pore and arises from the existence of 
metastable gas-like and/or liquid-like states. It is note- 
worthy that no calculations have been yet performed for 
cylindrical pores closed at one end. 

Macroscopic descriptions based on Cohan's model are 
currently used to describe the adsorption in non intercon- 
nected mesoporous media and in particular to determine 
the Pore Size Distribution (PSD) via the Barrett- Joyner- 
Halenda (BJH) method [§. 

However, the experimental validation of Cohan's 



model has never been checked since the synthesis of or- 
dered mesoporous media with non interconnected pores 
open at one or at both ends was unsuccesful for a long 
time. Using porous silicon, we report in this paper, for 
the first time, on an experimental investigation of Co- 
han's predictions. For this check, this mesoporous ma- 
terial exhibits two essential properties: (i) the tubular 
pores are non interconnected and do not exhibit narrow 
sections and (ii) it exists under two forms with pores 
closed at one end when the porous layer is supported 
by the substrate or open at both ends when the layer is 
removed from the substrate. 

Porous silicon is prepared by electrochemical etching of 
Si single crystal in a HF solution. The morphology of the 
porous Si layers depends on the type and concentration of 
the dopant atoms 0. The adsorption studies presented 
here have been performed with porous layers prepared 
from highly boron doped [100] Si substrate. The poros- 
ity of the layers is controlled by the current density and 
HF concentration, and the thickness by the duration of 
the anodic dissolution. The pore volume is determined 
with high precision by weighting the sample before and 
after the electrochemical dissolution. To calibrate the 
porosity and the thickness of the layer we need a third 
weighing after the dissolution of the porous part in a 
NaOH solution. It is noteworthy that the porous layers 
are perfectly reproduced. The porous layers under study 
had a porosity of 51% and the length of the pores is 21 
/im. 

The morphology of these layers (see Fig. |j](a) and 0(b)) 
can be described as honeycomb-like structure with pores 
perpendicular to the Si substrate separated from each 
other by Si single crystal walls (inset Fig. |l|(b)) of con- 
stant thickness (~ 5 nm). 

The pores, with sections of polygonal shape, are of the 
same length and closed at one end by the Si substrate. At 
the end of the porous layer formation which takes about 
20 minutes, an increase of the current density leads to 
the electropolishing regime during which the Si walls at 
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the bottom of the pores are dissolved so that the porous 
layer comes off the Si substrate. It is noteworthy that 
such a process takes only a few seconds and thus that 
the PSD of the porous layer is not changed. Playing 
with this possibility, porous Si layer or porous Si mem- 
brane can be obtained with pores open at one or at both 
ends, respectively. 
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FIG. 1. Bright field TEM images in plane view (a) and 
cross-section (b) of a porous Si layer with 51% porosity 
prepared from p + type (~ 3.10 -3 r2.cm) [100] Si substrate. 
Observation axis [100]. Pores (white) are separated by Si 
walls (black). The pore density is 3.2. 10 11 /cm 2 . Inset la: 
PSD. Inset lb: transmission electronic diffraction pattern. 



We have shown that there is no lateral interconnection 
between the pores in these layers. The method involves 
deposition by evaporation at ambient temperature, of an 
aluminum layer, 500 nm thick, on a part (A c and B c ) of 
the top of two porous samples A and B, following by a 
thermal annealing under neutral atmosphere of argon at 
450°C for half a hour to improve the contact between the 
Al deposit and the top of the layers. We have shown by 
Rutherford BackScattering analysis (RBS) with 2 MeV 
a particles that, the Al deposited on the samples is well 
localized on the top of the pores and constitutes a cap. 
One of the two samples (A) was thermally oxidized in 
2 enriched in 18 (99%), at 300°C, 12 mbar, for one 
hour. This treatment leads to the formation of a thin 
silicon oxide (~1 nm) on the silicon walls. The ls O con- 
tents in the porous layers under the Al cap (A c and B c ) 



and beside it (A Q ) were determined by the nuclear re- 
action ls O(p,a) 15 N at a proton energy near the narrow 
resonance at E p = 629 keV (r = 2 keV). For the mea- 
surements under the Al caps the energy of the protons 
beam energy was increased to take into account the en- 
ergy loss of the incident protons in the Al layer. The 
nuclear reaction spectra corresponding to A Q and A c are 
shown in Fig. |[ The a peak from A Q indicates of the 
presence in the porous layer of 18 O atoms coming from 
the thermal oxidation. Due to the presence of Al cap, 
the a peak from A c is expected at an energy lower by 
100 keV than that from the a peak from A Q . The counts 
from A c and B c which are one order of magnitude lower 
than that from A Q are equal, and hence corresponds only 
to a background coming from nuclear reactions with sili- 
con, boron atoms, etc. This experiment reveals that this 
medium is composed of non interconnected single pores, 
the top of which is in direct contact with the gas reser- 
voir during adsorption. To our knowledge, this is the 
only porous material for which the absence of intercon- 
nection has been experimentally proven. 





500 




400 






c 


300 


3 




a 




U 




200 




100 




2 2.5 3 3.5 4 

Energy (MeV) 

FIG. 2. a spectra from 18 0(p,a) 15 N: E p = 629 keV. 
White circles: porous Si layer after thermal oxidation in 
18 02 (A ). Black circles: part of the layer protected by Al 
cap (A c ). The arrow indicates the energy of the expected 
a peak from A c . 

Another important property of this medium is the con- 
stancy of the Si walls thickness revealed by the TEM 
plane view performed on a porous layer thinned down to 
1000 A, i.e. about 10 times the size of the pores. This 
is a strong argument against the presence of narrow sec- 
tions inside the pores, the presence of which would lead, 
through shadow effects, to an inhomogeneous apparent 
thickness of the Si walls. 

The PSD of the porous Si layers have been estimated 
from a numerical treatment of the TEM images. The 
negative TEM photographs are digitized by means of a 
CCD monochrome camera. Then, the value of the thresh- 
old to make the image binary is chosen to reproduce the 
porosity of the sample. The analysis of this binary image 
allows one to obtain the cross section area and perimeter 
of each pore. The PSD, which corresponds to cylindri- 
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cal pores having the same surface area as the polygonal 
pores, is shown in inset of Fig. [j](a). The mean diameter 
of the pores is 13nm with a standard deviation of ±6nm. 

Adsorption measurements were made by a volumet- 
ric technique using a commercial apparatus under sec- 
ondary vacuum of around 10~ 6 torr. Nitrogen adsorp- 
tion isotherms were measured up to the saturating vapor 
pressure of bulk nitrogen (P ) and at 77 K. Porous Si lay- 
ers were introduced in a cell just after their preparation. 
The samples were outgassed at ambient temperature for 
a few hours. This outgassing does not change the nearly 
complete passivation of the inner surface by a monolayer 
of SiHi^.3 groups which takes place during the electro- 
chemical etching and which is stable for a few days after 
the formation M. 
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where Vl is the molar liquid volume, a is the liquid-gas 
surface tension at temperature T, r is the radius of the 
pore and t(P) the thickness of the adsorbed layer at P. 
In the case of a cylindrical pore closed at one end, sur- 
face thermodynamics considerations predict a preferen- 
tial adsorption on the bottom of the pore; a hemispheri- 
cal meniscus is then present provided that there is perfect 
wetting, so that the condensation occurs at a pressure Ph 
given by the following modified Kelvin equation: 
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FIG. 3. (a) N2 adsorption isotherms at 77K for porous Si 
with pores open at one (black diamonds) or at both (white 
diamonds) ends, (b) Isotherms calculated by introducing 
in Eqs. (y) and (Q) the experimental t(P) values (see text) 
and the PSD shown in Fig. [j]. 

Figure |^(a) shows N 2 adsorption isotherms at 77 K cor- 
responding to porous Si layers with pores open at one or 
at both ends, the PSD of which is represented in inset of 
Fig. |l](a). These two curves exhibit wide and asymmet- 
rical hysteresis loops with shapes corresponding to type 
H2 with reference to the 1985 IUPAC classification (or 
type E in that of de Boer p). 

The adsorbed amount of gas sharply increases between 
0.55 and 0.8 Po and then reaches a plateau region where 
all the pores are filled with the dense phase. The steep 
desorption process occurs over a pressure range about 
three times narrower than that of the adsorption. 

On the assumption that the density of the dense phase 
equals that of the bulk liquid nitrogen, we found that 
the pore volume extracted from the value of the adsorbed 
amount on the plateau region is in a very good agreement 
with that measured by gravimetry. 

Cohan's model predicts that, for a cylindrical pore 
open at both ends, the condensation takes place at a 
pressure P c given by an equation similar to the classical 
Kelvin equation where the curvature of the meniscus is 
replaced by a cylindrical curvature : 



The evaporation of the wetting fluid proceeds through 
the same hemispherical meniscus in both cases i.e. at 
the same pressure Ph- 

Thus, according to Cohan's model : (i) the adsorption 
phenomenon should be reversible in pore open at one end 
while it should be irreversible in pore open at both ends; 
(ii) the adsorption branch for a pore open at both ends 
should be located at a pressure (P c ) higher than that 
for a pore closed at one end (Ph); (iii) for a given size 
distribution of cylindrical pores open at both ends, the 
adsorption branch should be steeper than that of the des- 
orption branch (provided that the diameter exceeds 3-4 
nm) and the shape of the hysteresis loop becomes rather 
symmetrical as the pore diameter increases, as shown by 
numerical calculations based on Eqs. ([j]) and (||) in which 
we have introduced a currently used t(P) proposed by 
Halsey |]ll| , |l2| or the t(P) reproducing our experimental 
values prior to the capillary condensation. 

To extract the latter t(P) values, we have prepared 
a porous Si layer with a PSD having a mean value as 
high as possible (about 50 nm) in order to shift the cap- 
illary condensation pressures to the highest possible val- 
ues. The obtained isotherm exhibits a hysteresis loop 
located above 0.8Po- We have thus determined the nu- 
merical values of t(P) up to O.8P0 by distributing homo- 
geneously at each pressure the adsorbed amount among 
the pores, the PSD of which was extracted as explained 
above. In order to draw the t(P) above O.8P0, we have 
chosen to extrapolate the polynomial law which fits the 
experimental t(P). 

The isotherm calculated by introducing in Eqs. ([!]) and 
(g) the t(P) and the PSD shown in inset Fig. |l](a) is rep- 
resented Fig. ||(b). 

We will successively discuss the predictions (i) and (iii) 
which arc in large discrepancy with the experimental re- 
sults (Fig. ||(a)); the prediction (ii) which is in agreement 
with experiment and finally, the large shift between the 
absolute position of the experimental isotherms and the 
calculated ones. 
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(i) The adsorption-desorption cycle observed in the 
case of pores closed at one end is irreversible. The in- 
terconnectivity related pore blocking effect described by 
Mason to explain the origin of the hysteresis phe- 
nomenon could not occur because the porous Si layers 
are composed of non interconnected pores. Moreover, 
we reject the 'ink bottle' effect, previously introduced by 
Everett fH}, as a possible explanation of the hysteresis 
loop we observe on the adsorption isotherms because of 
the absence of narrow sections inside the pores. 

(iii) The hysteresis loops observed in the porous layer 
and porous membrane exhibit the same shape of type H2. 
We have observed, by carrying out adsorption isotherms 
on porous Si layers with a PSD having mean values vary- 
ing from f 3 nm up to 50 nm, that the position on the pres- 
sure axis of the corresponding adsorption and desorption 
branches are shifted towards Po whereas the shape of the 
hysteresis loops remains of type H2. 

Such an asymmetrical shape is not expected by the 
single pore models. Indeed, numerical calculations based 
on Cohan's model or on DFT (l^] predict a symmetri- 
cal shape of the hysteresis loop for such large pores in 
contradiction with what we observe. Furthermore, al- 
though porous Si exhibits non-interconnected pores, the 
hysteresis shape (type H2) we observe is the same as 
those systematically obtained for highly interconnected 
mesoporous solids such as controlled pore glass or Vycor 
[15|. From a theoretical point of view, both macroscopic 
[ J and microscopic [|l3| approaches described this shape 
as a signature of the disorder introduced by the inter- 
connectivity between pores. We must conclude that a 
hysteresis loop of type H2 is not characteristic of inter- 
connected porous materials and we thus reject a previous 
conclusion made by Ball and Evans according to which 
a de Boer's shape of hysteresis is a consequence of inter- 
connected network effect. 

We believe that the desorption process we observe re- 
sults from another effect. At the end of the adsorption 
process, when all the pores are filled by the dense phase, 
the external surface of the porous layer is covered by 
a film which connects the pores. A possible interpre- 
tation of the observed delay to the emptying of pores 
could be the absence of formation of a concave menis- 
cus on top of each pore. Thus, question arises whether 
metastable states exist which would correspond to a flat 
interface between the film and gas. In such an idea, 
during the desorption process, the film thickness could 
decrease homogeneously until the system jumps to the 
stable state corresponding to a concave meniscus on top 
of each pore. The steepness of the desorption branch 
could be explained if the jump is sufficiently delayed un- 
til a pressure at which the emptying of the major part of 
pores occurs. This idea needs of course to be checked by 
numerical calculations. 

Note that the analysis of the desorption branch using 
BJH method yields a PSD with a width 0.5 nm about 



which is very far from that of the actual PSD (± 6nm). 
It is conventionally admitted jl6| that the desorption pro- 
cess in porous materials occurs at coexistence of the high 
and low dense phases. On such an assumption, a steep 
desorption branch is interpretated as the signature of the 
presence of pores of unique size ^J§. Our results obvi- 
ously show that it is not necessarily the case. 

(ii) The sign of the relative positions of the two adsorp- 
tion branches (Ph < P c ) is in agreement with Cohan's 
model. These experiments reveal two different scenari 
for the filling of pores open at one or at both ends. Pref- 
erential adsorption in bottom of closed pores is confirmed 
by analytical calculations performed in an oblique corner 
formed by the intersection of two planar surfaces [H . 
We are investigating pore closure effect on adsorption in 
cylindrical geometry by means of GCMC simulations and 
DFT calculations. 

As shown in Fig. ||(b), the curve which corresponds 
to the modified Kelvin equation (2) is largely outside of 
the experimental hysteresis loops. For pore sizes as large 
as in our samples, this equation should account for the 
equilibrium phase transition. Indeed, microscopic calcu- 
lations based on DFT performed by Ball and Evans [HJ 
have shown that, in the case of a distribution of cylin- 
drical pores open at both ends centered on 7.2nm, which 
is twice smaller than the mean size of our pores, the co- 
existence curve is located at a pressure which is close to 
that predicted by equation (2). To our knowledge, it is 
the first time that, for so large pore sizes, such a large 
discrepancy is found (f9) . This may be due to the use of 
a cylindrical geometry in the calculations which do not 
account for the actual polygonal shape of pores. GCMC 
simulations and calculations based on the DFT are in 
progress. 

In summary, we have studied adsorption phenomenon 
in porous Si medium, a well ordered material with tubu- 
lar pores without narrow sections and for which we 
have experimentally proven the absence of interconnec- 
tion. Hysteresis loops of type H2 are observed, and that, 
whether the pores are open at one or at both ends. These 
results are in contradiction with Cohan's model and show 
that the only observation of hysteresis loop of type H2 
is not the signature of the presence of interconnections 
in porous material as it is generally admitted. The steep 
desorption process cannot be described by the macro- 
scopic Kelvin equation in which we introduce the exper- 
imental values of the thickness of adsorbed amount prior 
to the capillary condensation. We note that the applica- 
tion of the widely used BJH method to extract the PSD 
from this branch is not systematically justified. 
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